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We consider the phase diagram of two-flavor quark matter under neutron star constraints for two 
nonlocal, covariant quark models within the mean field approximation. In the first case (Model I) 
the nonlocality arises from the regularization procedure, motivated by the instanton liquid model, 
whereas in the second one (Model II) a separable approximation of the one-gluon exchange inter¬ 
action is applied. We find that Model II predicts a larger quark mass gap and a chiral symmetry 
breaking (CSB) phase transition line which extends 15-20 % further into the phase diagram spanned 
by temperature (T) and chemical potential (p). The corresponding critical temperature at p = 0, 

T c (0) ~ 140 MeV, is in better accordance to recent lattice QCD results than the prediction of the 
standard local NJL model, which exceeds 200 MeV. For both Model I and Model II we have consid¬ 
ered various coupling strengths in the scalar diquark channel, showing that different low-temperature 
quark matter phases can occur at intermediate densities: a normal quark matter (NQM) phase, a 
two-flavor superconducting (2SC) quark matter phase and a mixed 2SC-NQM phase. Although in 
most cases there is also a gapless 2SC phase, this occurs in general in a small region at nonzero 
temperatures, thus its effect should be negligible for compact star applications. 

PACS numbers: 12.38.Mh, 24.85.+p, 26.60.+c, 97.60.-s 


I. INTRODUCTION 

Within the last decade the investigation of the thermodynamics and phase structure of strongly interacting matter 
has been driven by the results of the experimental programs with ultrarelativistic heavy ion beams at CERN-SPS and 
BNL-RHIC Q as well as by the unprecedented quality of data from lattice QCD simulations Q. A new picture of the 
state of matter created in these experiments has emerged, according to which the physical nature of the sought-for 
quark-gluon plasma (QGP) is aperfect liquid of strongly correlated hadron-like resonances rather than an ideal gas of 
quasifree quarks and gluons 0, 3 . The persistence of a strong, nonperturbative coupling in the QGP can be supposed 
as a prerequisite of a successful modelling of the QCD phase diagram. Turning to the domain of finite chemical 
potentials, guidance from lattice QCD is limited to p < T (where Taylor expansion techniques can be applied), and 
experimental programmes such as CBM at FAIR Darmstadt are yet in the stage of preparation with a planned begin 
of operation in 2014. Therefore, predictions for the QCD thermodynamics at low temperatures and high baryon 
densities p > T, where the critical point and the regions of color superconducting quark matter are expected in the 
QCD phase diagram [f|, have to be developed within effective models for nonperturbative QCD and to be tested 
against observational constraints from neutron stars [g]. For a recently developed testing scheme we refer to @]. 

After the early discussion on color superconducting dense quark matter 0 was revived for nonperturbatively strong 
couplings within effective quark models §j, a great variety of possible diquark pairing patterns and corresponding 
phases has been explored 3 and estimates for the order of magnitude of the pairing gaps have been given. One of 
the central questions for phenomenological applications in compact stars, where electric and charge neutrality has 
to be imposed EEL concerns the number of active flavors. It turns out that at low temperatures one obtains a 
sequential melting pattern of the light and strange quark chiral condensates, which is rather insensitive to the details 
of the four-momentum dependence of the interaction, but crucially dependent on whether also the strange quark 
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mass is determined sclfconsistently (see III for an early work within the covariant chiral quark model). The present 
’state-of-the-art’ for the description of color superconducting phases is represented by the fully selfconsistent mean 
field three-flavor Nambu-Jona-Lasinio (NJL) model calculations of Refs. 0 II 0 The results of these studies 
indicate that under compact star conditions, i.e. color and electric charge neutrality together with /?— equilibrium, 
the two-flavor color superconducting (2SC) phase is favored over the three-flavor color-flavor-locking (CFL) one. This 
outdates earlier results along the lines of, e.g., Ref. 0 where the strange quark mass has been set to a constant, 
small value. Even if the third quark flavor occurs at not too high densities to be in principle realized in compact star 
interiors, the star configurations with CFL quark cores turn out to be hydrodynamically unstable 0 . For a more 
detailed introduction to the recent status of dense quark matter in the NJL model, see Il7l . 

Studies of neutral 2SC quark matter consider also the presence of a so-called gapless (g2SC) phase El , which is 
found to occur at intermediate temperatures and chemical potentials. Going beyond the NJL theory, results arising 
from a noncovariant nonlocal quark model 00 with Gaussian regulator functions show that the 2SC phase is 
not present in asymmetric quark matter for standard values of the diquark coupling IHl . In this model, however, 
for strong diquark couplings one does find a 2SC phase, together with a g2SC and a mixed normal quark matter 
(NQM)-2SC phase in which the electric charge neutrality is satisfied only globally 000 . 

As the structure of the quark matter phase diagram within NJL-type models is sett ling^ it is important to investigate 
the situation in the case of effective models that go beyond those used in Refs. 0, G313 0> 0 j i n the sense that 
they include fully covariant nonlocal interactions. Nonlocality arises naturally in the context of several quite successful 
approaches to low energy quark dynamics as, for example, the instanton liquid model |25l and the Schwinger-Dyson 
resummation techniques ( 23 . The same happens in lattice QCD 1271. I t has been also argued that nonlocal models 
have several advantages over the local ones (e.g., the NJL model l28ll and its generalizations). Indeed, nonlocal 
interactions regularize the model in such a way that anomalies are preserved l29l and charges properly quantized, 
the effective interaction is finite to all orders in the loop expansion and there is no need to introduce extra cut-offs. 
Soft regulators such as Gaussian functions lead to small next-to-leading order corrections 30, 0], etc. This type of 
models has been successfully used to investigate meson 0000 and baryon 0,0, 0T properties at vanishing 
temperature and chemical potential. The phase diagram of isospin symmetric matter has also been studied within 
this context 000 0 , 0 ]. The aim of the present work is to extend these analyses to the case in which 
compact star conditions are imposed. We will consider two-flavor versions of the models, since (as mentioned below) 
effects arising from the strange quark can be safely neglected in the region of intermediate chemical potentials to be 
covered. The investigation of the question whether the results of the above referenced studies will remain qualitatively 
unchanged after the inclusion of correlations beyond the mean-field level goes beyond the scope of the present work 
and deserves a separate study. It should be noticed, however, that in the non-strange case considered here the relevant 
color superconducting condensate transforms as a singlet representation of the global SU(2)i ® SU(2)ij chiral group. 
Thus, no global symmetry is broken by the 2SC ground state, which implies that there are no Goldstone bosons that 
could further condense. 

The paper is organized as follows. In Section II we introduce the model and derive some useful analytical expressions. 
In Section III we present the numerical results obtained for the case of a Gaussian regulator, considering different 
ratios between the coupling constants. Then, in Section IV we discuss the features of the obtained phase diagrams, 
and in Section V we present our conclusions. 


II. FORMALISM 

Let us begin by stating the Euclidean action for the nonlocal chiral quark model in the case of two light flavors and 
anti-triplet diquark interactions, 

Se = j d 4 x ^(x) {-ifl + m) i>(x) - jj f M (x)j f M (x) - y \j a D {x)} ] j a D (x) j . (1) 

Here m is the current qua rk mass, which is assumed to be equal for u and d quarks. The nonlocality can be introduced 
now in different ways 0]. In what follows we will work within two alternative scenarios, that we call “Model I” and 
“Model II”, in which the mesonic current j’m(^) and the diquark current jo{x) in Eq. Q are given by nonlocal 
operators. In the case of Model I 0 , 0 ], the effective interactions are based on an instanton liquid picture of QCD. 
The nonlocal currents read 

j f M (x) = J d 4 y d 4 z r(y - x) r(x - z) fry) r f ${z) , 

3d( x ) = J d 4 y d 4 z r(y - x) r(x - z) $c{y ) * 75 ^ 2 A a ip(z ), 


( 2 ) 
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where we have defined ipc{x) = 7274 4> T {x) and T'f = (1, i^r), while t and A a , with a = 2,5, 7, stand for Pauli and 
Gell-Mann matrices acting on flavor and color spaces, respectively. 

On the other hand, Model II jd^Eh'j arises from a separable form of the effective one-gluon exchange (OGE) picture. 
In this case the nonlocal currents ]m,d(x) are given by 

3m( x ) = J d A z 9{z) ip(x + |) T f ip(x - '-), 

3d(x) = J d 4 z g(z) i> c (x + |) n 5 T 2 \ a ip(x - |) . (3) 

The functions r(x—y) and g{z) in Eqs. 0 and © are nonlocal regulators characterizing the corresponding interactions. 
It is convenient to translate them into momentum space. Since Lorentz invariance implies that they can only be 
functions of p 2 , we will use for the Fourier transforms of these regulators the forms r(p 2 ) and g{p 2 ) from now on. 

The current-current interaction in Eq. 0 is quite common to effective theories for low-energy QCD such as the 
NJL model. The momentum dependence introduced by the functions r(p 2 ) and g{p 2 ) is a generalization of that 
model, and can be chosen in such a way that the dynamical mass function of the quark propagator reproduces some 
features obtained in lattice QCD analyses. It could entail, e.g., the possibility of complex conjugate mass poles of 
the quark propagator, in agreement with results arising from the Schwinger-Dyson equation approach to QCD. The 
index structure of the vertices Tf and also the ratio of the coupling constants H/G can be obtained from a Fierz 
rearrangement of e.g. the OGE interaction, as detailed in na (see also references therein). For the OGE, as well as 
for the instanton liquid model, one obtains H/G = 0.75. However, since a precise derivation of the effective couplings 
from QCD is not known, there is a significant theoretical uncertainty in this value. In fact, so far there is no strong 
phenomenological constraint on H/G , except for the fact that values larger that H/G ~ 1 are quite unlikely to be 
realized in QCD, since they might lead to color symmetry breaking in the vacuum. We will leave this ratio as a free 
parameter, analyzing the results obtained for values lying within a range from 0.5 to 1. 

The partition function of the system at temperature T and quark chemical potentials p,f c is given by 

Z = j V//V1I) e ~ s ^fc,T) ) ( 4 ) 

where the Euclidean action is obtained from Eq. 0 by going to momentum space and performing the replacements 


P4 > Cd n 5 



(5) 


Here p^ is the fourth component of the (Euclidean) momentum of a quark carrying flavor / and color c, and ui n are 
the Matsubara frequencies corresponding to fermionic modes, ui n = (2 n + 1 )ttT. We are assuming here that quark 
interactions depend on the temperature and chemical potentials only through the arguments of the regulators. Note 
that, as required for the applications we are interested in, we have introduced in Eq. 0 different chemical potentials 
for each quark flavor and color. 

To proceed it is convenient to perform a standard bosonization of the theory. This procedure is described, e.g., 
in Refs. 0, 0|. Thus, we introduce the bosonic fields <7, n a and A a , and integrate out the quark fields. In what 
follows, we work within the mean field approximation (MFA), in which these bosonic fields are replaced by their 
vacuum expectation values 7 r a = 0, a and A a . Moreover, we adopt the usual 2SC ansatz A 5 = A 7 = 0, A 2 = A. We 
have assumed other possible condensates to be negligible, following previous analyses carried out in the NJL model 
framework mu Within the MFA, and employing the Nambu-Gorkov formalism, the thermodynamical potential 
per unit volume can be written as 


n MFA 


T 

~ In Z MFA 


, ]A£ 
2 G 2 H 



In det 


T 


S-\a, A) 


( 6 ) 


Here the inverse propagator S 1 (cr, A) is a 48 x 48 matrix in Dirac, flavor, color and Nambu-Gorkov spaces, given by 
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where we have used the definitions 


pfc = ( u>n T * Pic , P ) , 

(8) 

E/c = rn + a h(pj c ,pj c ) , 

(9) 

A = A h(pi r ,p dr ) , 

(10) 


with / = u,d and c = r, g, b. The functions h(p, q) have been introduced in order to have a common notation for both 
Model I and Model II. One has 


h(s, t ) 


r(s 2 ) r(t 2 ) 



(Model I) 
(Model II) 


( 11 ) 


We have taken into account that, as we will see below, the usual 2SC ansatz implies /if r = Hfg- In Eq. 0, entries 
with subindices ur and dr are intended to be multiplied by an 12x2 matrix in rg space, while r| stands for a T 2 Pauli 
matrix acting in this space. 

The determinant of S' -1 can be analytically evaluated, leading to 


n MFA = 


=.2 


2G 2 H 


- T E 


d 3 p \ A c\ 2 


n =—oo 


(2tt) 3 ^ T 8 ’ 


where 


A c — (pi c + E 2 C ) (p dc + A, dc 2 ^ + (1 — 6bc ) A 2 [A 2 + 2 ( pi c ■ p dc ) + 2 E uc E^ c ] . 


( 12 ) 


(13) 


For finite values of the current quark mass, Q MFA turns out to be divergent. The regularization procedure used here 
amounts to define 


up, = n MFA - n free + n { r e 


+ (reg) > 


(14) 


where O free is obtained from Eq. HU by setting A = cr = 0, and is the regularized expression for the thermo¬ 

dynamical potential of a free fermion gas, 


^free = _ 2 
a ^(reg) 


T f Win [l + e -(^ +m2 -^)/ T l + In [l + e -(V^+™ 2 +M/c) /Tl j _ 

d 1 n ' f,c 1 1 11 11 


(15) 


The mean field values a and A are obtained from the coupled gap equations 

dVLf FA , 

(reg) = A (1 - 16 HT D ud ) = 0 , 

= a - 4 G T (S ud + S du ) = 0 , 

LIU 

where we have defined 

OO 

Dij — ddji — Re ^ ) 


dA 

r l QMFA 

da 


d 3 p 


d (Puri P dr ) 


A 2 + (pf r ■ p jr ) + E ir E* 


Ar 


S%j = Re ^ 


n =—oo 


■ d 3 p 
(27r) 3 


+ , S ir (p. r +E* r 2 ) + A 2 E* r E, b 

2 h(pt r ,p+) -1-j-- + h(p+,p+) 


pi +£U 


(16) 

(17) 

(18) 
(19) 


So far we have introduced different chemical potentials for each quark flavor and color. However, not all of them 
are in general independent quantities. For the description of quark matter in the core of neutron stars, we require 
the system to be color and electric ch arg e neutral (for a further discussion on the issue of color neutrality and 
color singletness we refer to Refs. 0 |4flj). Thus, within the previously introduced 2SC ansatz, only one color- 
dependent chemical potential is needed 1471 . and the g,f c can be written in terms of only three independent quantities: 
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the baryonic chemical potential hb 7 the quark electric chemical potential HQ q and the color chemical potential h g. 
Defining h = hb/3, the corresponding relations read 

2 1 

Hur — Hug — H + T^HQq + —^8 

1 1 

Hdr — Hdg ~ H ~ T^HQq + 3M8 

2 2 

Hub = H + r^HQq - gM8 

1 2 

Hdb = H- ^HQq ~ 3M8 ( 20 ) 

Now, in the core of neutron stars, in addition to quark matter we have electrons. Thus, within the mean field 
approximation for the quark matter, and considering the electrons as a free Dirac gas, the full grand canonical 
potential is given by 


n fuU 


= n 


MFA 

( re g) 


<r 


( 21 ) 


where 


^ = - 1^ ’ (22) 

fi e being the electron chemical potential. For simplicity we have neglected here the electron mass. 

In addition, it is necessary to take into account that quark matter has to be in j3— equilibrium with electrons 
through the beta— decay reaction 


d —> u + e + v e . 


(23) 


Thus, assuming that antineutrinos escape from the stellar core, we must have 


l^dc l^uc — MQq — Me • 


(24) 


If we impose the requirements of electric and color charge neutrality, /x e and H8 become fixed by the conditions of 
vanishing electric and color densities: 


PQtot 




(pfr + Pfg 2pfb) ~ 0 , 

f=u,d 


(25) 


where 


on _ dn e 

dHe dfJLe 


On d^eg) 

d Hfc d/ifc 


(26) 


Consequently, in the physical situation we are interested in, for each value of T and h we should find the values of A, 
a, He and hs that solve Eqs. m and G3> supplemented by Eqs. (Ehfll and l|25|> . 


The electron density can be evaluated analytically, leading to 


On the other hand, from Eqs. ClJ, JI3 and m the quark densities can be expressed as 


(27) 
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(28) 


Then a straightforward calculation leads to 
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with 


dA c 

d^uc 


2 ( iu) n + Hue) (iPdc A- ^-“dc ) (1 "h 2 E UC S UC ) 

— 2 (1 She) A' 2 £^ c S uc ( i u) n -\- Hue) A i ^n Hdc 

dHt’Pde) 


- 4* (1 - 5 bc ) A A 


A" + (Pte-Pde) + 


dtA 


t=pl 


(30) 


where we have defined EV = a dh(t,t)/dt 2 | 2= + 2 . 

The corresponding expressions for pdc are obtained by simply exchanging u and d and taking the complex conjugate, 
while the expressions for can be easily obtained from by setting a = A = 0. Finally, p(reg)/c given by 


free _ 9 
P(reg)/c - 1 


d 3 p 

( 2^)3 


1 + exp 


\Jp 2 + m 2 - Hfe 
T 


1 + exp 


yf- 


Hfe 


T 



(31) 


III. NUMERICAL RESULTS 


In this section we present our numerical results, showing the features of the phase diagram and the behavior of 
relevant physical quantities for Models I and II. According to previous analyses carried out within nonlocal scenar¬ 
ios |42|, the results are not expected to show a strong qualitative dependence on the shape of the regulator. Thus we 
will concentrate here on simple and well-behaved Gaussian regulator functions, taking (in momentum space) 

r{p 2 ) = exp(-p 2 /2A 2 ) (Model I) (32) 

g(p 2 ) = exp(— p 2 /A 2 ) (Model II) (33) 

Here A is a free model parameter, playing the role of an ultraviolet cut-off. We have chosen a different normalization 
for Models I and II in view of the relation between the respective regulating functions [see Eq. CEDI, which determine 
the low T and h phenomenology. 


A. Parameterization 

For definiteness, for both Models I and II we choose here input parameters m, A and G which allow to reproduce the 
empirical values for the pion mass m„ = 139 MeV and decay constant /„ = 92.4 MeV, and lead to a phenomenologically 
acceptable value for the chiral condensates at vanishing T and Hfe- For Gaussian regulators, taking into account the 
chosen normalization of the cut-offs, it is seen that within the MFA both models lead to the same expressions for 
the considered physical quantities at T = h = 0. However, this is not the case when one goes beyond the MFA [isj . 
In particular, the expressions for the pion mass and decay constant are different (they are still coincident only in 
the chiral limit), and it is necessary to use different sets of input parameters. The parameters considered here for 
Model I are m = 5.14 MeV, A = 971 MeV and GA 2 = 15.41, while for Model II we have taken m = 5.12 MeV, 
A = 827 MeV and GA 2 = 18.78. With these sets we get for both models a phenomenologically reasonable value for 
the chiral condensate, namely (OlgglO) 1 / 3 = —250 MeV. The remaining free parameter is the coupling strength H in 
the scalar diquark channel. In order to fix its value by hadron phenomenology at zero T and h , one would have to 
solve the Faddeev-type equations for baryons as three-quark bound states which result from the quantization of chiral 
quark models of the type considered in the present paper after bosonization in meson and diquark channels (see fTil 
and Refs, therein; for more elaborate recent calculations of nucleon properties see [33, ESI ) ■ Although in principle 
this is possible, we refrain from fixing H by hadron phenomenology within the present exploratory study of the quark 
matter phase diagram and rather choose different values for the coupling ratio H/G in the range from 0.5 to 1. 


B. Order parameters and phase transitions 

For fixed values of the temperature T and the chemical potential h (h = Hb/3), the mean field values d and A, as 
well as the chemical potentials He and hs , can be numerically obtained from the gap equations m and im together 
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FIG. 1: Behavior of the mean fields a and A and the chemical potentials fi e and ps for Models I (left) and II (right) as a 
function of the chemical potential, for three different values of the temperature. The curves correspond to the case H/G = 3/4. 
Full lines correspond to T = 0, dashed lines to T = 40 MeV and dotted lines to T = 100 MeV. In the case of T = 0, lines 
marked with stars and dots correspond to the 2SC and NQM phases respectively. 


with the conditions of f3— equilibrium and vanishing color and electric charge densities, Eqs. EH) and EH- Let us 
begin by considering the case H/G = 0.75, which is motivated by various effective models of quark-quark interactions. 
Our results for ft, A, fi e and /is are shown in Fig. Q where we plot these quantities as functions of fi for different 
representative values of the temperature. Left and right panels correspond to Models I and II, respectively. 

For T = 0 (solid lines), at low chemical potentials the system is for both Models I and II in a chiral symmetry 
broken phase (CSB), where quarks acquire large dynamical masses. By increasing the chemical potential one reaches 
a first order phase transition, in which the chiral symmetry is approximately restored, and a certain volume fraction of 
the quark matter undergoes a transition to the 2SC phase coexisting with the remaining normal quark matter (NQM) 
phase. The chemical potential fi e (which for T = 0 vanishes in the CSB region) also shows a discontinuity across the 
transition. The new 2SC-NQM mixed phase is a way in which the system realizes the constraint of electric neutrality 
globally: the coexisting phases have opposite electric charges which neutralize each other, at a common equilibrium 
pressure. In its simplest realization, this mixed phase is treated within an approximation in which Coulomb and 
surface energies are neglected (see Ref. |m|). For color superconducting quark matter this realization of charge 


























neutrality has been considered e.g. in Ref. [24[ for the NJL model and in Ref. [2lJ for the instantaneous nonlocal 
quark model. The discussion of inhomogeneous mixed phases, which are not yet fully understood, crucially depends 
on the assumptions for the surface tension and charge screening effects (see e.g. Ref. |52l )t and goes beyond the scope 
of the present investigation. On the other hand, following Refs. 23j and |53f], we have imposed color neutrality 
as a local constraint. This is based on the fact that the color Debye screening length is expected to be short and 
comparable to the inter-particle distance in the regime of interest. As a consequence, fis turns out to be different in 
the two components of the mixed phase. However, it should be kept in mind that this chemical potential is in fact 
an effective quantity that has to be introduced in these kind of models in order to account for the effect of the gauge 
fields. Namely, as argued in Refs. HEIM superconducting quark matter is expected to be automatically color 
neutral in QCD. As expected, the growth of the color chemical potential fig in the 2SC component of the mixed phase 
is approximately proportional to that of the corresponding A, which governs the amount of breakdown of the color 
symmetry due to quark pairing. 

When the temperature is increased (see dashed curves in Fig.^ corresponding to T = 40 MeV), the mixed phase 
is no longer favored and the system goes into a pure 2SC phase. For T = 40 MeV, this shows up as a second order 
transition in the case of Model I, and a first order transition in the case of Model II. Now, for both models, when 
one moves along the first order transition line from T = 0 towards higher temperatures, one arrives at a triple point 
(3P). At this point the CSB and 2SC phases coexist with a third NQM phase, in which the chiral symmetry is 
approximately restored and there is no color superconductivity. Finally, if T is still increased, one reaches an “end 
point” (EP) where the first order transition from CSB to NQM phases becomes a smooth crossover. The behavior 
of the dynamical masses and the electric chemical potential fi e along this smooth transition is shown in Fig. ^ see 
curves corresponding to T = 40 MeV (Model I) and T = 100 MeV (dotted lines, Models I and II). 

C. Quark matter phase diagrams 

The described features of the phase diagrams for Models I and II can be visualized in the graphs shown in Fig. [3 
where we plot the transition curves on T — fi diagrams for different ratios H /G, and show the regions corresponding 
to the different phases and the position of triple and end points. In the graphs, solid and dotted lines correspond 
to the mentioned first order and crossover transitions, respectively (in the case of the crossover, the transition point 
can be defined by considering the maximum of the chiral susceptibility Uli). Between NQM and 2SC regions we 
find that in all cases there is a second order phase transition, which corresponds to the dashed lines in the diagrams 
of Fig. 0 Close to this phase border, the dashed-dotted lines in the graphs delimit a band that corresponds to the 
so-called gapless 2SC (g2SC) phase. Here, in addition to the two gapless modes corresponding to the unpaired blue 
quarks, the presence of flavor asymmetric chemical potentials fide — l^uc ^ 0 gives rise to another two gapless fermionic 
quasiparticles 0 . Although the corresponding relations cannot be derived analytically owing to the nonlocality of 
the interactions, the border of the g2SC region can be numerically found. This is done by determining whether for 
some value of \p\ the imaginary part of some of the poles of the Euclidean quark propagator vanish in the complex 
P 4 plane. From the graphs it is seen that this g2SC band may become relatively significant for low values of H/G. 
However, it never extends up to zero temperatures, therefore this should not represent a robust feature for compact 
star applications. 


IV. DISCUSSION 

Let us discuss some qualitative features of the curves displayed in Figs. ^ and |3 On one hand, for both models 
the 2SC phase region becomes larger when the ratio H/G is increased. This is not surprising, since H is the effective 
coupling governing the quark-quark interaction that gives rise to the pairing. In any case, as a general conclusion 
it can be stated that, provided the ratio H/G is not too low, these nonlocal schemes favor the existence of color 
superconducting phases at low temperatures and moderate chemical potentials. Indeed, for the parameters considered 
here, there is no 2SC phase only in the case of Model II, H/G = 0.5. This is in contrast with the situation in e.g. the 
NJL model 00, where the existence of such a phase turns out to be rather dependent on the input parameters. In 
addition, our results are qualitatively different from those obtained in the case of the noncovariant nonlocal models [2lj . 
where above the chiral phase transition the NQM phase is preferable for values of the coupling ratio H/G < 0.75. In 
those models, a color superconducting quark matter phase can be found only for H/G ~ 1 . 

It is also interesting to compare our results with those obtained for isospin symmetric quark matter. For the same 
parameter sets, the corresponding phase diagrams for H/G = 0.75 are shown in Fig. [3] By comparing them with those 
of Fig. [3 it can be seen that the 2SC region becomes reduced when one imposes color and electric charge neutrality 
conditions. This is indeed what one would expect, since the condition of electric charge neutrality leads in general 
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FIG. 2: Phase diagrams for Models I (left) and II (right) for different values of H/G. Full and dashed lines indicate first and 
second order phase transition curves respectively, dotted lines correspond to crossover-like transitions, and dashed lines delimit 
the gapless 2SC band. Different phases are denoted as NQM (normal quark matter phase), CSB (chiral symmetry broken 
phase) and 2SC (two-flavor superconducting phase), while the regions marked as “Mixed” correspond to the NQM-2SC mixed 
phase. Note that in Model I, for H/G = 1 there is a narrow region in which one has 2SC while the system is still in a chiral 
symmetry broken phase. EP and 3P denote the end points and triple points respectively. 
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to unequal u and d quark densities, disfavoring the u-d pairing. We notice, however, that the effect is relatively 
small, and the positions of triple and end points as well as the shape of the critical lines remain approximately un¬ 
changed. Concerning the shape of the chiral phase transition line Tcsb(m)j one observes at intermediate temperatures 
50 MeV < Tcsb 100 MeV, which are relevant for the future CBM experiment, an approximately linear behavior. 
This can be seen as an interpolation between a convex shape obtained in NJL or bag models and a concave shape for 
confining Dyson-Schwinger equation models, see e.g. Ref. M- It is remarkable that thus in the nonlocal covariant 
models presented here a similarity with confining quark models occurs and that the chiral/deconfinement transition 
line in the phase diagram resembles very closely the positions of freeze-out parameters in heavy ion collisions. 
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FIG. 3: Phase diagrams for symmetric matter corresponding to Models I (left) and II (right). Here we have taken H/G = 0.75 


Finally, we have studied the dependence of the phase transitions on the model parameters, changing the input value 
of the chiral condensate within a phenomenologically reasonable range 220 MeV < —(0|^<j r |0) 1,/3 < 280 MeV. From this 
analysis, it is seen that the qualitative features of the phase diagrams are not significantly modified. In particular, it 
is seen that one finds in general color superconducting phases at low temperatures and moderate chemical potentials, 
for intermediate values of the ratio H/G. In addition, the values for the critical temperature at fi = 0 are quite stable, 
yielding about 120 MeV for Model I and 140 MeV for Model II. This would favor the description given by Model II, 
in which the result is closer to the values provided by lattice calculations. 


V. CONCLUSIONS 

We have considered the phase diagram of two-flavor quark matter under neutron star constraints for two nonlocal, 
covariant quark models within the mean field approximation. In the first one the nonlocality was due to the regu¬ 
larization procedure, motivated by the instanton liquid model (Model I), whereas in the second model a separable 
approximation of the one-gluon exchange interaction was applied (Model II). Although for the Gaussian regulators 
considered in this work the vacuum gap equations are identical, both models differ in their fluctuation spectrum and 
therefore in their parameters (current quark mass m, coupling strength G, UV cutoff A), which have been fixed by 
using input values for the pion mass, the pion decay constant and the chiral condensate. As result of the numerical 
evaluation of the corresponding gap equations at finite temperature and chemical potential, we have obtained that 
Model II predicts a larger quark mass gap and a chiral symmetry breaking phase transition line which extends 15-20 % 
further into the T — fj, plane when compared to Model I. The prediction for the critical temperature at /i = 0 in Model 
II, Tcsb ^140 MeV, is closer to the results of recent lattice QCD simulations than the prediction of both Model I 
and the well-known local NJL model. Considering different values for the coupling strength in the scalar diquark 
channel, we have found that under neutron star constraints different low-temperature quark matter phases can occur 
at intermediate densities: normal quark matter (NQM), pure superconducting (2SC) quark matter and mixed 2SC- 
NQM phases. The critical temperature for the 2SC phase transition is a rising function of /i in the case of Model 
I whereas it is rather independent of ^ for Model II, due to the different /x dependences associated with the scalar 
diquark gaps. A band of gapless 2SC (g2SC) appears at the border between 2SC and normal quark matter. At large 
values of the H/G ratio it is given by a tiny strip in the phase diagram, and grows broader at low diquark couplings. 
However, the g2SC region does not reach zero temperatures, thus it should not represent a robust feature for compact 
star applications. Our present investigation has been limited to the mean field approximation and the neglect of the 
strange quark flavor. None of these approximations is expected to be crucial for compact star applications, and the 
detailed study of their impact on our results has been left for further development of the present approach. 
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